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FINAL REPORT (W911NF0810061): A Combined 
Experimental/Computational Investigation of Nanoelectromechanical 
Switches - Understanding and Eliminating Prevalent Failure Modes 

ABSTRACT 

The International Technology Roadmap for Semiconductors (ITRS 1 ) identifies emerging technologies with 
the potential to sustain Moore's Law. A necessary succession from planar CMOS to non-planar/dual-gate 
CMOS, and ultimately to novel architectures such as carbon nanotube- (CNT)-based nanoelectromechanical 
systems (NEMS) is envisioned. The ITRS also identifies critical roadblocks currently precluding advances 
beyond CMOS. Primary among the roadblocks to NEMS are poor reliability and manufacturing challenges. 
Here we investigate, both computationally and experimentally, the prevalent failure modes of CNT-based 
NEMS that hamper reliability. We first identify their point of onset within the design space, highlighting the 
extremely limited region in which failure is currently avoided. We use dynamic multiphysics models to 
elucidate the underlying causes of failure, and then show that the usable design space expands dramatically 
when employing novel electrode materials such as diamond-like carbon. We then demonstrate the efficacy of 
this solution through numerous successive actuation cycles without failure and applications to volatile 
memory operations. Finally, we develop a probe-based nanomanufacturing scheme by which to scale up 
manufacturing of these robust devices. Ultimately, these advances will be critical in delivering a broad class 
of robust, battle-ready NEMS sensors and electronics. 


1. OBJECTIVES 

NEMS technologies are widely pursued for various applications where technological demands require reduced 
component size and cost, while dramatically increasing device performance and functionality. These competing 
factors present significant challenges in device design. Within this project, we develop a platform CNT-NEMS 
architecture, focusing on understanding the failure modes that are prevalent to this class of devices and 
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ultimately on overcoming them to allow robust application to memory systems and sensors. Upon successfully 


overcoming the prevalent failure modes, we turn to develop manufacturing schemes that will allow construction 
of large-scale arrays of these robust devices. Specifically, our objectives are to: 

1. Investigate emerging CNT-NEMS technologies, which are likely to become the next generation of 
electronic devices and sensors. 

2. Perform a dynamic analysis of the device and investigate its performance in the context of applications 
such as NEMS switches, sensors, random-access memory elements, and logic devices. 

3. Experimentally identify prevalent device failure modes and their point of onset within the design space. 

4. Employ multiphysics modeling to study the devices' dynamic pull-in response and explain the underlying 
mechanisms for experimentally-observed failure modes. 

5. Computationally predict conditions (device geometry and pull-in voltage) under which device failure is 
likely and design new architectures devoid of these failure modes. 

6. Develop manufacturing schemes to fabricate large-scale arrays of these robust devices. 

7. Train graduate students and post-docs in a multidisciplinary environment. 

Ultimately, accomplishing these objectives will demonstrate the feasibility of transitioning CNT-NEMS from the 
lab bench to the battle field in the form of truly robust, highly functional devices. 

2. INTRODUCTION & SIGNIFICANCE 

The immense potential of CNT-based NEMS is emergent in theoretical and experimental demonstrations of 
up to 100-GHz switching 2 , low leakage and high ON-OFF ratios 3 , and outstanding current-carrying capacity 4 ' 5 . To 
date however, individual demonstrations of performance such as these have been a primary focus, with limited 
reports of repeated actuation beyond a few cycles 2 " 3,6 " 7 . This is due in large part to critical challenges in device 
reliability and manufacturing. As a consequence of their ultra-high performance, these devices experience 
extreme (and coupled) mechanical, electrical, and thermal loads during normal operation. Furthermore, 
previously insignificant forces (e.g., van der Waals interactions) begin to dominate at this scale. These factors 
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combine to create a number of failure modes that are prevalent among reported devices. Thus finding a solution 


to suppress or eliminate the onset of failure is critical to the widespread integration of CNT-based NEMS and will 
represent a significant milestone along the ITRS path. 

Given solutions to these failure modes, scaling up manufacturing to produce large-scale arrays of CNT-NEMS 
poses another significant ITRS roadblock. The serial methods used to produce the aforementioned devices are 
in general not conducive to scaling up. For example, electron beam lithography is perhaps the most commonly 
used method to define electrical contacts to randomly dispersed carbon nanotubes. However, this inherently 
serial process requires involved imaging, alignment, and pattern generation steps tailored specifically to each 
individual nanostructure. To date, no method exists to create large-scale arrays of well-ordered individual 
carbon nanotubes that is readily integrated with other microfabrication processes to construct functional 
devices. 

In this work, an electrostatically-actuated CNT switch 8 ' 13 serves as a platform to study prevalent failure 
modes in CNT-based NEMS and methods by which to scale up manufacturing. This switch is chosen because it 
shares operating principles (and thus failure modes and manufacturing challenges) with many reported devices, 
making the findings of this study broadly applicable 1415 . The device consists of a multi-walled CNT of length L 
that is fixed at one end and cantilevered a distance H above an electrode (Figure la). The electrical circuit is 
completed by a voltage source, V, and an external resistor, R f . When the applied bias exceeds a critical "pull-in 
voltage" (V > l / P i), the CNT accelerates toward the electrode, closing the switch. Subsequent reduction of the 
applied bias results in re-opening of the switch 8 . The electrical domain can be represented by an equivalent 
lumped-element circuit (Figure lb) comprised of the voltage source and external resistor in series with, in 
parallel, the parasitic capacitance of the system C 0 , the capacitance of the CNT cantilever relative to the 
electrode C C nt, and the tunneling/contact resistance, R C nt, when the CNT approaches and makes contact with 
the electrode. Note that C CW rand /? C /vrwill depend upon H (see Section 4.1). The total current is denoted as / tota i, 
while the current through the CNT is / CNT . U is the potential across the CNT, where U = V when / tota | = 0. 
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Figure 1 Overview of cantilevered CNT switch used as a platform to study failure modes, (a) Schematic of switch consisting of a CNT 
cantilevered over a metal thin film electrode, (b) Equivalent lumped-element circuit for the device, (c) Nanomanipulaton-based 
experimental setup used to recreate the device geometry shown in (a), enabling real-time incremental variation of device geometry and 
simultaneous electromechanical characterization in situ the SEM. The voltage source V is the only component exterior to the SEM 
chamber, (d) Scanning electron micrograph showing a CNT mounted on the tip of a nanomanipulator probe approaching the electrode. 
Scale bar is 1 pm. 


3. EXPERIMENTAL IDENTIFICATION OF FAILURE MODES 

A parametric study was first conducted to identify prevalent failure modes and their point of onset within 
the device design space. To test a large number of device geometries spanning the design space in a more 
statistically-significant study, an in situ nanomanipulation-based technique was adopted over conventional 
microfabrication processes to construct the devices. A 3-axis nanomanipulator was used to position 
cantilevered CNTs over an electrode (Figure lc-d, see Section 7.1), effectively recreating the platform device 
(Figure la). This allowed real-time prescription of device geometry with simultaneous voltage application (l/), 
measurement of total current (/ tota |), and imaging of device operation and failure via SEM. For example, a device 
of a given length (!) could be repeatedly characterized with incrementally decreasing gap (/-/) by using the 
manipulator to progressively step the CNT closer to the electrode. In this initial parametric study, a gold-coated 
substrate was used as an electrode as reported devices ubiquitously use metal thin film electrodes ' ' ' . 
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Two common modes of failure were observed in the initial experiments (Figure 2). For relatively long 


cantilevers and small CNT-electrode gaps, irreversible stiction preventing re-opening of the switch was observed 
(Figure 2a-c,f). The corresponding measured l to tarV response (Figure 2a) shows a well-defined jump in current 
upon pull-in, followed by a linear decrease to zero as the applied voltage is ramped back down. This linear 
behavior is characteristic of maintained Ohmic contact between the CNT and electrode. Imaging the device 
following actuation confirmed that the CNT indeed remained stuck to the electrode after the applied bias was 
fully released (Figure 2c). This irreversible stiction has been widely reported and is known to be the result of 
strong van der Waals interactions between the CNT and electrode 2 ' 6 ' 16,19 ' 20 . When these adhesive interactions 
exceed the elastic restoring force in the deformed cantilever, the switch does not re-open, even when the bias is 
fully removed. In fact, nanotube-substrate interactions have been shown in some cases to be sufficiently strong 
to radially collapse the CNT 21 ' 22 . In cases where the device is operated in atmosphere, liquid adhesion adds 
additional adhesive energy to be overcome. 
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Figure 2 Characteristic failure modes, (a-c) Failure by irreversible CNT-electrode stiction. (a) Characteristic measured / toto r'/ response for 
a device exhibiting irreversible stiction upon pull-in. Voltage was applied in a single 0-15 V triangle wave. (b,c) Scanning electron 
micrographs showing the CNT cantilever before and after the actuation cycle respectively. In (c), the CNT clearly remains stuck to the 
electrode after removing the applied bias. Scale bars are 250 nm. (d-e) Failure by incremental loss of CNT length, (d) Characteristic 
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measured / toto/ -V response for a device that does not suffer from irreversible stiction. Voltage was applied in a single 0-35 V triangle 
wave. As the voltage is ramped down after pull-in, the current falls slightly non-linearly until pull-out. Inset shows detail of the pull-out 
event, after which the current is slightly negative due to gradual discharging of the capacitances, C 0 and C CNT , as the applied voltage 
continues to decrease, (e) Plot of CNT length (normalized by length before first actuation cycle) following 10 successive actuation cycles, 
(e) Length of CNT cantilever after successive actuation cycles exhibiting incremental loss. Length is normalized by the initial length before 
the first cycle, (f) Map of device geometries tested in the L-H design space and the corresponding failure modes observed. Lines indicate 
analytically-derived boundaries between failure modes (see Section 4.4). Below a critical elastic restoring force, irreversible stiction is 
expected as adhesive forces dominate. Above a critical current density (p = p c ), CNT ablation is expected upon pull in due to Joule 
heating. 

For relatively short cantilevers with larger CNT-electrode gaps, incremental shortening of the CNT cantilever 
was observed with successive actuation cycles (Figure 2d-f). Here the characteristic l to tarV curve exhibits a sharp 
drop in current (inset of Figure 2d) as the voltage is ramped down after pull in, signifying successful re-opening 
of the switch. However, inspection of the CNT after successive actuation cycles revealed an incremental loss of 
length (Figure 2e). Similar incremental losses have been reported, as well as more extreme cases of complete 
loss of the CNT 8,23 ' 25 . Two possible causes are envisioned. During the highly-dynamic pull-in event, the CNT 
experiences large mechanical stresses as it impacts the electrode at high velocity, which could fracture the CNT. 
Alternatively, Joule heating resulting from large currents could ablate the nanostructure, resulting in similar loss 
of length 23 ' 26 . 

4. DYNAMIC MULTIPHYSICS SIMULATIONS TO EXPLAIN FAILURE 

To investigate the underlying cause for the observed incremental loss and test potential solutions, a 
dynamic multiphysics finite element model of the device was developed using ABAQUS/Explicit and appropriate 
user-defined subroutines. This model, which includes consideration for electrostatic forces and van der Waals 
interactions between the CNT and substrate, enabled detailed analysis of the time history of the mechanical 
stresses and electrical current in the CNT. 

4.1 Development of Model 
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In the model, the sum of the attractive van der Waals force and electrostatic force on the CNT were treated 
as a body force and implemented through the VLOAD subroutine. The repulsive van der Waals force was used 
to simulate the interaction between the CNT and the substrate and implemented through the VUINTER 
subroutine. 

The dynamic response of the CNT is described as 27 ' 28 


P A J^ + c ^ + EI lhc* ~ Clelec + ^ vdw ' ( 1 ) 

where p = 2.27 g/cm 3 is the density of the CNT, A is the cross-section area of the CNT, c is the damping constant 
of the system, E is the elastic modulus and / is the moment of inertia. On the right-hand-side of (1), q eiec and 
Rvdw represent the load per unit length on the CNT due to electrostatic and van der Waals interactions, 
respectively. The damping constant of the system is given as 


c — 


nApco 0 

Q ’ 


( 2 ) 


where Q is the quality factor and u) 0 is the natural frequency of the CNT, which can be described as 


29 


_ P 2 I El 

0)0 l 2 Ja p ' 


(3) 


where (3 = 1.875. The quality factor is taken as 250, which is estimated from the vibration of a CNT cantilever 
driven by alternating electrical field 29 . 

The electrostatic interaction between a CNT and an electrode has been studied in previous work 12,30 . The 
electrostatic force applied per unit length of the CNT is given as 

Qeiec = n£ 0 U , ^ + ^ ( 4 ) 

lr(r+2R)-a cosh 2 ^l +—j 


where £ 0 is the permittivity, V is the applied potential on the CNT, r is the gap between the CNT and the 
substrate, R is the external radius of the CNT, and f c is a function which considers the concentrated charge 
effect on the free end of the CNT cantilever 11 


fc = 0.85 [(H + «) 2 «] 1 / 3 5(x - x tip ) (5) 

in which H is the initial gap between the CNT and substrate. 
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The van der Waals interaction between the CNT and substrate is calculated based on the Lennard-Jones 
potential 31 . For the interaction between two atoms with distance, r , the Lennard-Jones potential is 

*(r) = £\-% (6) 

in which £12 and C 6 are material properties. For carbon-carbon interaction 32 , C 12 = 24.1 keV A 12 and C 6 = 15.2 
keV A 6 . The potential between an atom and a semi-infinite solid body (with atom density, p x ) can be obtained 
from the integration 

= C C 1°oo 3 Pirfxdydz (7) 

where D is distance between the atom and the semi-infinite solid. Derivative of the potential with respect to D 
yields the force acting on the atom as 

/«» = £ = ¥{£-&}■ < 8 > 

Thus, the force acting on a unit volume of atoms (with atom density, p 2 ) is 

= 5-zB}- <=» 

In the model, and p 2 are set as the atomic density of graphene, 114/nm 3 . 

Therefore, the total attractive force between the CNT and the bottom electrode is the sum of the attractive 
van der Waals interaction (the first term on the right-hand-side of (9)) and the electrostatic force, and can be 
written as 


q+ (D) = eim^ + . 

^ + v J 5 2D 4 


—TT£qV Z 


: (1 + fc) 


2 Rt Jr(r+2R)-acosh 2 (l+^j 

The repulsive van der Waals interaction (the second term on the right-hand-side of (9)) is used to define the 
contact between the CNT and the bottom electrode. Therefore, the normal contact pressure on the CNT is given 


where pf represents the number of atoms per unit area. For graphene, pf = 38/nm 2 . 
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Since the bias charge on the CNT is time dependent, it is necessary to investigate the equivalent lumped 


circuit as shown in Figure lb to obtain the charge on the CNT in real time. The equivalent circuit yields 

R f (_C 0 + C c „ T )% + (R/-^ + 1 ^ ;+ l)u = V ( 12 ) 

in which Rf is the feedback resistance, C 0 is the capacitance between the two electrodes, R C nt ' s the tunneling 
resistance of the CNT to the electrode. In vacuum, the tunneling resistance between the CNT and bottom 
electrode can be described as R C nt = R c ex p(A5), where R c is the contact resistance, A is a constant 22.8 nm' 1 
and 6 is the gap between the CNT tip and the substrate 33 . In the subroutine, the bias on the CNT, U, is updated 
at each time step by solving (12) numerically. 

4.2 Description of Simulations 

To investigate the cause of the observed incremental shortening of the CNT with repeated actuation and 
test possible solutions, the model described above was setup as follows. In all reported simulations, the CNT 
cantilever had a length of 1.5 pm, an outer radius of 25 nm, and a wall thickness of 25 shells. The assumption of 
25 shells is based on the reported average number of shells for multi-walled CNTs obtained from the same 
source 34 . The CNT was modeled as an orthotropic linear elastic material 35 . Four-node shell elements were used 
with an element size of 2 nm to investigate the buckling of the CNT during pull-in. The gap between the 
centerline of the CNT and the electrode was 100 nm. A voltage was applied in a step function of 15 V beginning 
at t = 0 sec to actuate the device. 

For cases with a gold electrode, the electrode was modeled as an elastic/plastic material with a Young's 
modulus of 78 GPa and a Poisson's ratio of 0.44. A low yield stress of 10 MPa was chosen to count for the 
softening of Au due to Joule heating. For cases with a diamond-like carbon electrode (tetrahedral amorphous 
carbon (ta-C), see Section 5), the electrode was modeled as a purely elastic material with a Young's modulus of 
87 GPa and a Poisson's ratio of 0.22. For both cases, a value of /? f = 1 GQ was used which, in the case of the gold 
electrode, is representative of the external resistor placed in the circuit (experimentally-measured to be 992 
MQ). In the case of ta-C, the geometry of the ta-C electrode itself was designed to have a resistance of 1 GQ 
between the point of CNT contact and its connection to the voltage source (see Section 7.2). Thus no physical 
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resistor was needed in the circuit, with the resistance of the ta-C electrode/trace taking its place instead. An 


experimentally-measured value of C 0 = 50 pF was used for the parasitic capacitance (see Methods). 
Experimentally-measured values of contact resistance were used for both cases as well, with R c = 50 kQ for gold 
and R c = 0.6 GQ for ta-C (see Section 7.3). To isolate the effect of the electrode on the device behavior, all other 
parameters were held constant. Eight-node linear brick elements with reduced integration were used for 
modeling the bottom electrode. 

4.3 Results of Simulations 

In the simulations, a step input is given at the voltage source, causing the CNT to pull in. As the CNT strikes 
the electrode (Figure 3a), a stress wave propagates from the point of contact toward the constrained end, and 
reflects back (Figure 3c-f). Additionally, minor buckling is observed momentarily along the CNT (Figure 3e). 
Though significant, the peak stress observed is approximately one order of magnitude less than the reported 
failure stress for CNTs obtained from the same source 27 . Thus mechanical fracture alone is an unlikely cause of 
the observed incremental shortening with repeated actuation (Figure 2e). 



6.1 ns 


6.3 ns 


'0 


Figure 3 Results of dynamic multiphysics finite element simulation investigating the cause of incremental CNT loss, (a) CNT tip- 
electrode gap (H) versus time in response to a step in voltage applied at t = 0. The response is shown for devices with gold (solid blue 
line) and ta-C (dashed red line) electrodes. Insets show snapshots of deformed CNT profile as it deflects toward the electrode, (b) 
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Current through CNT (/ CNT , see Figure lb) versus time for devices with gold and ta-C electrodes. Numbered points along the curves 


correspond in time to those in (a), (c-f) Snapshots of stress in tip of the CNT during CNT-electrode impact. 

Instead of catastrophic mechanical stresses, the dynamic multiphysics analysis reveals a transient spike in 
current through the CNT (Figure 3b) as it accelerates toward the electrode. This spike, which approaches the mA 
range, is the result of rapid charge dissipation from the capacitances, C 0 and C CN t (Figure lb), to the substrate. 
The corresponding current density (see Section 4.4.1) is more than sufficient to ablate the CNT 4 . In fact, similar 
observations of incremental loss have been reported in CNT field emitters operated at similar current 
densities 23 ' 25 . Plotting a line of constant current density in the L-H design space equivalent to the critical current 
density required to burn a multi-walled CNT with 25 shells (Figure 2f, see Section 4.4.1) separates well the 
experimentally-tested cases in which incremental loss was and was not observed, reaffirming the predictions of 
the model. 

During the initial discharge, / CN t » Itotai (Figure 3b). Once the stored charges have been fully depleted and 
the switch is fully closed, the system reaches a significantly lower steady-state current flow (/ tota i = Icm = 
V/(R C NT+fy), see Figure lb). Due to the computational cost of the dynamic model, it was not possible to simulate 
a sufficiently long time period to reach this steady state. However, this steady state current has been previously 
predicted analytically for similar devices 8 and is measured in the present experiments (/ tota i. Figure 2a and Figure 
5a). Here it is important to note that, based upon the equivalent lumped-element circuit (Figure lb), the spike in 
/ CNT is not expected to appear in the measured l tota rV response (Figure 2d). Because the spike arises due to 
discharge of the capacitances, it remains within the /? c/VT C 0 loop and does not contribute directly to the 
measured / tota |. 

4.4 Predicting the Onset of Failure 

To test the suspected modes of failure as predicted by the FEM model against the experimental results and 
thus validate its predictive power, analytically-derived boundaries between the onset of the various failure 
modes were plotted over the experimentally-observed cases (see Figure 2f and Figure 5e). These boundaries 
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were derived based upon the assumed modes of failure predicted by the model. Thus their match with the 
experimental results will support the validity of the model. 

4.4.1 Incremental loss ofCNT length due to Joule heating-induced ablation 

The results of the multiphysics FEM modeling suggest that a significant discharge current occurring when 
the CNT touches the gold electrode is sufficient to ablate the CNT. This is the result of low CNT-electrode 
contact resistance leading to a fast time constant for dissipation of charge stored in the system. The reported 
threshold current density ( p c , also known as breakdown current density) is on the order of 10 8 A/cm 2 for CNTs 
and graphene 5,36 . Similar observations of incremental loss from CNTs used as field emitters have been reported 
for high current densities 23 . 

When the CNT contacts the bottom electrode, the total charge stored in the system will dissipate rapidly 
due to the low contact resistance. Given that \J~V PI at the time of contact, the maximum current density on the 
CNT can be approximated by 


Pi = 


Vpi 


(14) 


Rcnt 2nR ext hN' 

where R ext is the external radius of the CNT, h is taken as 0.34 nm (the interlayer spacing between CNT shells), 
and N is the number of shells in the multi-walled CNT. The pull-in voltage is expressed as 12 


V PI = 


0.85 


1+ 


8 H 2 


1 + 


2.55 [Rext(H + RextP] ^ 



(15) 


where L is the tube length. Combining (14) and (15) and setting p t = p c yields the condition when the outer 
shell tends to burn when the CNT contacts the Au electrode. 

4.4.2 Irreversible stiction due to van der Waals interactions 



Figure 4 Schematic of stiction of CNT on bottom electrode due to van der Waals interaction. 
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As the length of the CNT cantilever increases, the CNT tends to stick irreversibly to the bottom electrode 


after contact. This is because the restoring spring force in the deformed CNT (which decreases with length) 
cannot overcome the van der Waals interaction with the substrate. Here the stiction length (the length of the 
CNT that is in contact with the substrate, see Figure 4) is denoted as L c . When in contact with the substrate, the 
CNT can thus be treated as a cantilever with a length of (L — L c ) whose left end (x = 0) is clamped and at the 
right end (x = L — L c ), the deflection is H with no rotation. The restoring spring force of the cantilever is 
expressed as 


12 EIH 
0 i-LcY' 


(16) 


The van der Waals interaction can be expressed as 


F V dW — a ^C' (17) 

where a is the van der Waals force per unit stiction length. Irreversible stiction occurs when F vdw > F re , yield 
a stiction condition of 


L c (L-L c ) 3 >^, (18) 

with the constraint that 0 < L c < L. The right-hand-side of (18) has maximum of 3 3 L 4 /4 4 at L c = L/4. 

Therefore, the van der Waals stiction yields a profile of constant critical restoring force in the L-H design 
space (Figures 2f and 4c in the manuscript) which follows H or L 4 . The coefficient in this quartic relation will 
depend upon the CNT-electrode material interaction, a. 


4.5 Using the Model to Redesign the Device Architecture 

Looking at the highly limited remaining region of the design space in which no failure is predicted (white 

region in Figure 2f), it is not surprising to find such limited reports of robust, repeated cycling of electrostatically- 

actuated CNT-based NEMS. Perhaps the simplest way to avoid ablation would be to reduce the CNT-electrode 

gap to decrease the required pull-in voltage (and thus the stored charge). However, this promotes stiction. To 

avoid this, other devices have employed more complex three-terminal electrode configurations in which a gate 

electrode at high bias is used to push or pull the CNT into contact with a drain electrode at relatively low bias 16 " 
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In this way the CNT makes contact only with the drain electrode that is at a relatively small potential 


17 , 37-39 

difference. However, this requires more complex accessing/multiplexing schemes to scale up and reduces the 
potential integration density. As such, reports of devices operating without failure or degradation through 
numerous cycles remain elusive with these architectures. 

Examining the equivalent circuit more closely (Figure lb), increasing R CNT would moderate the detrimental 
current spike by increasing the time constant for the charge dissipation. However, the contact resistance 
between CNTs and gold or other conventional metallic electrodes is generally low (measured to be in the kQ 
range in our experiments), resulting in rapid discharge and a large current spike at contact 26 . Unfortunately, 
simply increasing the external resistor R f would not have a similar effect as it is outside the RcntCo loop and thus 
does not affect the time constant of the discharge. 

Diamond-like carbon (tetrahedral amorphous carbon, ta-C) was next considered as an alternative electrode 
material. Here the contact resistance between the CNT and electrode was measured to be 0.6 GQ (see Section 
7.3), 5 orders of magnitude greater than for gold. Based upon the previous analysis, this should slow the charge 
dissipation upon pull-in, thus mitigating the current spike. Before implementing experimentally, this hypothesis 
was tested using the dynamic multiphysics finite element model. A device of the same geometry was re¬ 
simulated with an increased value of contact resistance (/? CNT = 0.6 GQ, see Section 7.3). The resulting mechanics 
of the switch closing are nearly identical, resulting in similar stresses. Figure 3a shows the CNT tip-electrode gap 
as a function of time in which the cases using gold and ta-C electrodes follow nearly identical paths. By contrast, 
the spike in current through the CNT drops dramatically as expected in the case of the ta-C electrode. Figure 3b 
compares the current profile through the CNT using the gold and ta-C electrodes. Here we see that the 
magnitude of the current spike is less than 2.5 pA (as compared to > 300 pA for gold), dropping the resulting 
current density well below the critical value for burning. As the CNT accelerates toward the electrode (point #3 
in Figure 3a-b), a small spike in / C n T is still observed due to the rapidly increasing capacitance with decreasing gap 
before contact which results in charges being pumped into the CNT. However, because of the increased time 
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constant for capacitance discharge, the stored charge is dissipated over a significantly longer time, resulting in a 


peak /cnt current closer to the steady-state measured current. 




o 50 100 

Time (Actuation Cycle) 



Figure 5 Suppressed failure modes and demonstrated device functionality using diamond-like carbon electrodes, (a) Characteristic / totar 
V curve showing well-defined ON/OFF behavior with significantly less stiction than devices with gold electrodes. Inset shows detail of the 
pull-out event, after which the current is slightly negative due to gradual discharging of the capacitances, C 0 and C CNT , as the applied 
voltage continues to decrease, (b) Current profile for 100 successive actuation cycles driven by a 0-35 V triangle wave at the voltage 
source. Inset shows a detail of cycles 46-50. The numbered data points correspond to the numbered positions in the / totar V curve in (a). 
(c,d) Demonstration of basic volatile memory operation, (c) Schematic of voltage input (l/ !N ) used to write '0' or T and ouput (1/ 0 ut)- (d) 
Voltage input and measured output for a series of read and write operations, (e) Map of devices tested in the L-H geometric design space 
with ta-C electrodes. Incremental shortening/ablation is eliminated, while the onset of stiction is greatly suppressed as compared to 
those with gold electrodes. For comparison, the dashed area represents the robust, failure-free region determined for devices using gold 
electrodes (Figure 2f). 


5. ELIMINATING FAILURE MODES THROUGH INFORMED DEVICE DESIGN 

Given these positive predictions, devices with ta-C electrodes were then tested using the same in situ 
experimental technique (see Section 7.1). The characteristic / to tai-V curves exhibited clean pull-in response and 
well-defined pull-out (Figure 5a) as compared to the gold electrodes (Figure 2a). The observation of sharper pull¬ 
out occurring earlier relative to the pull-in voltage suggests decreased CNT-electrode adhesion. This is reflected 
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in the map of failure modes (Figure 5c) where we see greatly suppressed onset of irreversible stiction along the L 


axis for a given H. This can be attributed to the relatively low surface energy and large contact angle of ta-C as 
compared to metals such as gold 40 ' 42 . 

Perhaps more important though is the complete elimination of any observed incremental loss of CNT length 
during repeated actuation when using ta-C electrodes (Figure 5e). Even in devices in which the CNT-electrode 
gap approached the total length of the CNT, no shortening was observed. This elimination of observed 
shortening supports the model prediction of significant reduction of the current spike. Rather than being bound 
by charge dissipation, the design space is now limited by a less severe geometric constraint in which the CNT- 
electrode gap (/-/) simply cannot exceed the length (/.) of the CNT. 

To highlight the efficacy of this solution, a single device was tested through 100 successive actuation cycles 
by applying a triangle wave at the input voltage. Imaging the device before and after the 100 cycles showed no 
observable loss of length. Throughout the series of actuation cycles, the measured current versus time (Figure 
5b) profile corresponds well to the current versus voltage measured before the cycles (Figure 5a). Given that 
the predicted mechanical stresses for the ta-C-based devices are similar to those of the Au-based devices, this 
successful high-cycle demonstration supports the model's hypothesis that mechanical stresses are not the cause 
of incremental CNT loss when Au electrodes are employed. It further supports the premise that charge 
dissipation-induced Joule heating is instead the key mechanism and demonstrates the efficacy of increasing 
CNT-electrode contact resistance in eliminating this failure mode. 

To further demonstrate the functionality of the devices given their improved robustness, basic memory 
operation was tested. Treating the voltage source (Figure lb and Figure 5c) as the input and the potential of the 
electrode as the output (Figure 5c) allows volatile storage of '0' and T states. Figure 5d shows a train of inputs 
and the corresponding measured output. To write 'O', 1/| N <V'pull-out causes the switch to open and l/ 0UT to go low. 
To write '1', l/| N > 1/ PU ll-in causes the switch to close and 1/ 0 ut to go high. To read or store, 1/pull-out<^n<^pull-in 
holds the switch in its current state. 
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6. SCALING UP MANUFACTURING 


To take full advantage of the robust CNT-based NEMS devices developed above, multiple devices operating 
in parallel are necessary. For example, simple OR logic requires two such two-state devices 43 . Beyond individual 
logic operations, the ability to make numerous devices in parallel will allow for highly multiplexed sensors and 
other electronic devices. As with commonly-used electron beam lithography processes, the in situ 
nanomanipulation technique used herein to investigate failure modes of CNT-NEMS is not conducive to scaling 
up. Instead, scaling up requires the ability to create large, well-ordered arrays of individual CNTs in parallel. 

¥ 


Figure 6 Schematic of nanoscale delivery using the Nano Fountain Probe 44 , (a) Liquid ink is stored in an on-schip reservoir and fed 
through enclosed microchannels to aptertured writing tips by capillary action. The inset shows a scanning electron microscopy (SEM) 
image of an aptertured tip (scale bar is 2 pm), (b) For direct-write nanopatterning, the tip is brough into contact with the substrate 
where an ink meniscus forms. 

To address this challenge within the objectives of this project, a unique probe-based nanomanufacturing 
tool called the Nano Fountain Probe (NFP 45 " 46 , Figure 6) was leveraged to create precise arrays of catalyst from 
which carbon nanotubes could subsequently be grown 47 . The NFP, an atomic force microscopy (AFM) tool 
developed in the laboratory of the PI, functions like a highly miniaturized fountain pen to deliver liquid 
molecular "inks" with sub-100-nanometer resolution. These inks are loaded into an on-chip reservoir and fed 
through a series of enclosed microchannels to apertured writing tips (Figure 6). When the tip is brought into 
contact with the substrate, the ink forms a stable meniscus, allowing transport to the substrate for deposition 
(Figure 6b). By taking advantage of the precise spatial resolution of the AFM by which the NFP positioning is 
controlled during patterning, extremely precise arrays can be formed. Previously-demonstrated applications 
include nanopatterning of gold 48 and diamond 44 nanoparticles, biomolecules 49 ' 50 , drugs 44 , and thiol 51 solutions 
on a variety of substrates. 
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Figure 7 Nano Fountain Probe-based atalyst patterning for large-scale carbon nanotube arrays, (a) Feature size (dot diameter and 
height) as a function of square-root dwell time, (b) Dot array of catalyst patterned on a silicon dioxide substrate using the Nano Fountain 
Probe. Scale bar is 4 pm. (c) Two dots of catalyst from which carbon nanotubes have grown. Scale bar is 1 pm. 

Arrays of catalyst were patterned on a variety of substrates using the NFP 47 , including silicon, silicon dioxide, 
silicon nitride, quartz, and platinum. The ability to pattern precise dot arrays (Figure 7b) of three different 
catalyst types was demonstrated (iron(lll) chloride, cobalt acetate, and ferritin). Feature size (dot diameter and 
height) showed strong dependence on the dwell time (Figure 7a), enabling precise control over the amount of 
catalyst deposited. Carbon nanotubes were then grown from the catalyst dots (Figure 7c) by chemical vapor 
deposition in a methane atmosphere at temperatures of 750 to 850 C. 



Figure 8 Schematic of a CNT-NEMS device fabricated using the Nano Fountain Probe-based catalyst patterning technique. First, all 
necessary electrodes, traces, vias, etc are constructed using conventional microfabrication techniques. As a final step, the Nano Fountain 
Probe is used to place the catalyst in desired locations and the carbon nanotubes grown in place. 

This unique ability to create well-ordered arrays is now being integrated with a complete process flow to 
construct large numbers of the robust devices developed in Sections 3-5. Within this process, conventional 
microfabrication techniques are first used to create all necessary electrodes, traces, vias, etc. The NFP is then 
used to place catalyst in the locations where carbon nanotubes are desired (Figure 8). As a final step, the chips 
are placed in the CVD furnace and suspended nanotubes grown in place from the catalyst. This method of 
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device fabrication is unique in that it not only allows parallel fabrication of numerous devices, but also minimizes 


the exposure of the carbon nanotubes to harsh processing steps commonly used in microfabrication (e.g., strong 
etchants which can introduce defects in the nanotube structure) by placing the nanotube growth step at the end 
of the process flow. In contrast, most e-beam lithography processes to create freestanding nanotube structures 
make use of carbon nanotubes randomly dispersed on a sacrificial layer. After defining electrical contacts, the 
sacrificial layer is removed using strong chemical etchants (e.g., hydrofluoric acid) to release the nanotubes. 
When combined with the findings of the previous sections to eliminate prevalent failure modes (i.e., through the 
use of diamond-like carbon electrodes), these results will be a significant advance toward battle-ready NEMS. 


7. SUMMARY AND SIGNIFICANT OUTCOMES 

In summary, the limited reports of successful repeated actuation of CNT-based NEMS using conventional 
metal electrodes are not surprising given the prevalence of irreversible stiction and ablation failure modes 
throughout the design space as identified in this work. A dynamic multiphysics model of the device revealed 
that, in addition to the well-understood irreversible stiction, another dominant failure mode occurs due to Joule 
heating-induced ablation. This is the result of rapid dissipation of stored charges upon closing of the CNT switch. 
Together, these failure modes drastically confine the useful region of the geometric design space. We then 
showed that the robust, failure-free region of the design space is dramatically expanded through the use of 
alternative electrode materials such as diamond-like carbon. Using these diamond-like carbon electrodes, we 
see not only reduced stiction, but complete elimination of the ablation failure mode. This enables an 
unprecedented number of actuation cycles and reliable function as volatile memory, providing motivation for 
further implementation of alternative materials in improving performance, robustness, and reliability of 
electrostatically-actuated NEMS. Finally, given the efficacy of using diamond-like carbon electrodes in 
eliminating failure modes, we developed a probe-based nanomanufacturing scheme to create well-ordered 
arrays of individual carbon nanotubes. This will be integrated with a microfabrication process flow to construct 
large-scale arrays of robust CNT-based NEMS to enable parallel operation of these switches for logic operations. 
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To conclude, the significant outcomes from this project are: 


1. Comprehensive identification of prevalent failure modes of CNT-based NEMS and their point of 
onset within the geometric design space using in situ electromechanical characterization. 

2. Development of a dynamic multiphysics model which accurately captures the electromechanical 
response of the devices. This model, which was validated against experimental results, reveals the 
underlying mechanisms for the experimentally observed failure modes while enabling facile testing 
of potential solutions to these failure modes. 

3. Dramatic improvement in device reliability by replacing commonly-used metal thin film electrodes 
with diamond-like carbon electrodes. This enabled hundreds of actuation cycles without failure and 
greatly expanded the usable region of the design space. This redesign of the device was based upon 
findings from the dynamic multiphysics model. A provisional patent application has been filed for 
this technology. 

4. Development of a probe-based scheme for patterning arrays of catalyst for subsequent growth of 
carbon nanotubes in large, well-ordered arrays. 

Ultimately, accomplishing these outcomes demonstrate the feasibility of transitioning CNT-NEMS from the lab 
bench to the battle field in the form of truly robust, highly functional devices. 
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8. METHODS 


8.1 In Situ Device Characterization 

In situ electromechanical device characterization was conducted in a Nova NanoSEM (FEI) scanning electron 
microscope (SEM) for simultaneous high resolution imaging and current-voltage measurements. Arc-discharge- 
grown multi-walled CNTs (n-Tec, Norway) were first dispersed in dry powder form on a copper TEM grid 
(Structure Probe, Inc). The grid was mounted on a Teflon block and placed in the SEM. SEM imaging was used 
to identify protruding CNTs of the desired length and diameter for device characterization (for all experiments, 
CNTs of consistent radius (21 ± 3 nm) were selected to allow direct comparison of results between devices of 
varying length and gap). A 3-axis piezoelectric nanomanipulator (Klocke Nanotechnik, Germany) mounted inside 
the SEM was then used to bring a sharpened tungsten probe (Micromanipulator Co., USA) into contact with the 
free end of the selected CNT. Electron beam-induced deposition of carbon was used to fix the CNT in place on 
the tip of the probe. This deposition technique has been shown to produce near-ideal Ohmic contacts to CNTs 52 . 
The nanomanipulator was then used to extract the cantilevered CNT and bring it into close proximity to the 
electrode (also mounted on the same Teflon block to electrically isolate it from the SEM stage). This effectively 
recreated the platform device geometry (Figure la) with the CNT-electrode gap being precisely controlled by the 
position of the nanomanipulator probe relative to the electrode (Figure lc). 

The probe and electrode were electrically addressable via a feedthrough in the SEM chamber wall. The 
circuit was completed by insertion of a 1 GO resistor (Ohmite Mini-Mox, 992 MQ actual measured resistance) in 
line between the electrode and SEM chamber feedthrough. An electrical bias was applied and current 
simultaneously measured using a parameter analyzer (Keithley 4200-SCS Semiconductor Characterization 
System). During l to ta\-V measurements, the electron beam of the SEM was blanked as it typically generated noise 
in the measurements as high as 0.5 nA. 

The CNTs were imaged before and after actuation to inspect for failure. CNT length and diameter, and CNT- 
electrode gap were measured from SEM images using ImageJ software (NIH). The CNT diameter was taken from 
the full width at half max of the intensity profile across the cross section of the CNT in the image. 
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8.2 Electrode Fabrication 


Gold electrodes were fabricated by depositing a 100-nm film of gold (with a 10 nm chromium adhesion 
layer) on a 200-nm-thick silicon nitride-coated silicon wafer by thermal evaporation. A 1-cm 2 piece was then 
cleaved from the wafer and fixed to the Teflon block as described above. Ta-C electrodes were fabricated by 
depositing a 140-nm film of ta-C on similar silicon nitride-coated silicon wafer. To define the proper total 
electrode resistance and eliminate the need for the external 1 GO resistor, long thin traces of ta-C were 
patterned in the ta-C film. First, a 70-nm aluminum film (with a 10-nm titanium adhesion layer) was deposited 
by evaporation over the ta-C and patterned by photolithography and liftoff. This was used as an etch mask to 
define the ta-C electrode shape. The exposed ta-C was etched through to the silicon nitride by reactive ion 
etching (RIE) using CF 4 /0 2 . The aluminum etch mask was then stripped by RIE using BCI 3 /CI 2 /He to re-expose the 
ta-C electrodes. Finally the ta-C sample was mounted on a Teflon block in place of the gold electrode. 

8.3 Measuring Parameters for Multiphysics Model 

The multiphysics models require knowledge of the parasitic capacitance and CNT-electrode contact 
resistances. These were measured experimentally to improve the quantitative accuracy of the models. To 
measure the parasitic capacitance, the nanomanipulator was first used to bring a CNT into close proximity with 
the electrode as before. The Keithley 4200-SCS system was then used to source a constant current of 30 fA 
while measuring the corresponding rise in voltage required to source this current. Assuming the parasitic 
capacitance is significantly larger than the CNT-electrode capacitance (C 0 » C CNT , where C CNT is on the order of 
10' 17 F 12 ) the parasitic capacitance was estimated based upon the time required to ramp from 0 to 1 V at 
constant current (C 0 = i-ht/hV). This yielded a parasitic capacitance of about 50 pF for both the Au- and ta-C- 
electrode systems. 

The measured parasitic capacitance was validated by comparison to the capacitance predicted based upon 
the magnitude of the slightly negative current observed after pull-out (insets of Figure 2d and Figure 5a). After 
pull-out, the switch is open again at which point the only current expected would be due to discharging of the 
parasitic capacitance in response to the continued decrease in applied voltage. Given a fixed value of 
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capacitance, this current is proportional to the rate which the voltage was ramped down (/' = C 0 -/W//\t). In these 


experiments, the voltage was ramped down at a rate of -2.2 V/s which, given a capacitance of 50 pF, would yield 
a current of -0.11 nA as the voltage is ramped down. This compares favorably to the experimentally measured 
currents of approximately -0.1 nA after pull-out, adding confidence to the values used for parasitic capacitance. 

The contact resistance was estimated first by taking the slope of the l to tarV curve after pull-in to calculate 
the total circuit resistance. From this, the total circuit resistance when the tungsten probe was pressed into 
direct contact with the electrode (no CNT present) was subtracted to yield the CNT-electrode contact resistance. 
This calculation is based upon the assumption that the probe-electrode contact resistance is significantly less 
(because the area of contact is orders of magnitude greater) than the CNT-electrode contact resistance. This 
yielded a CNT-electrode contact resistance of R c = 5-50 kQ for gold, and R c = 0.6 GQ for ta-C. Note that in the 
equivalent lumped element circuit (Figure lb), the combined tunneling/contact resistance term (/? CNT = R c 
exp(A5), where 5 is the gap between the CNT and substrate) approaches R c as 5 -> 0 (see Section 4.1). 
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